Pan-Cancer Analysis of lncRNA Regulation Supports Their Targeting of Cancer Genes in Each Tumor Context

Long noncoding RNAs (lncRNAs) are commonly dys-regulated in tumors, but only a handful are known toplay pathophysiological roles in cancer. We inferredlncRNAs that dysregulate cancer pathways, onco-genes, and tumor suppressors (cancer genes) bymodeling their effects on the activity of transcriptionfactors, RNA-binding proteins, and microRNAs in5,185 TCGA tumors and 1,019 ENCODE assays.Our predictions included hundreds of candidateonco- and tumor-suppressor lncRNAs (cancerlncRNAs) whose somatic alterations account for thedysregulation of dozens of cancer genes and path-ways in each of 14 tumor contexts. To demonstrateproof of concept, we showed that perturbations tar-geting OIP5-AS1 (an inferred tumor suppressor) andTUG1 and WT1-AS (inferred onco-lncRNAs) dysre-gulated cancer genes and altered proliferation ofbreast and gynecologic cancer cells. Our analysis in-dicates that, although most lncRNAs are dysregu-lated in a tumor-specific manner, some, includingOIP5-AS1, TUG1, NEAT1, MEG3, and TSIX, synergis-tically dysregulate cancer pathways in multiple tumorcontexts

Pretreatment Serum Concentrations of 25-Hydroxyvitamin D and Breast Cancer Prognostic Characteristics: A Case-Control and a Case-Series Study

Results from epidemiologic studies on the relationship between vitamin D and breast cancer risk are inconclusive. It is possible that vitamin D may be effective in reducing risk only of specific subtypes due to disease heterogeneity.In case-control and case-series analyses, we examined serum concentrations of 25-hydroxyvitamin D (25OHD) in relation to breast cancer prognostic characteristics, including histologic grade, estrogen receptor (ER), and molecular subtypes defined by ER, progesterone receptor (PR) and HER2, among 579 women with incident breast cancer and 574 controls matched on age and time of blood draw enrolled in the Roswell Park Cancer Institute from 2003 to 2008. We found that breast cancer cases had significantly lower 25OHD concentrations than controls (adjusted mean, 22.8 versus 26.2 ng/mL, p<0.001). Among premenopausal women, 25OHD concentrations were lower in those with high- versus low-grade tumors, and ER negative versus ER positive tumors (p≤0.03). Levels were lowest among women with triple-negative cancer (17.5 ng/mL), significantly different from those with luminal A cancer (24.5 ng/mL, p = 0.002). In case-control analyses, premenopausal women with 25OHD concentrations above the median had significantly lower odds of having triple-negative cancer (OR = 0.21, 95% CI = 0.08-0.53) than those with levels below the median; and every 10 ng/mL increase in serum 25OHD concentrations was associated with a 64% lower odds of having triple-negative cancer (OR = 0.36, 95% CI = 0.22-0.56). The differential associations by tumor subtypes among premenopausal women were confirmed in case-series analyses.In our analyses, higher serum levels of 25OHD were associated with reduced risk of breast cancer, with associations strongest for high grade, ER negative or triple negative cancers in premenopausal women. With further confirmation in large prospective studies, these findings could warrant vitamin D supplementation for reducing breast cancer risk, particularly those with poor prognostic characteristics among premenopausal women

Genomic, Pathway Network, and Immunologic Features Distinguishing Squamous Carcinomas

This integrated, multiplatform PanCancer Atlas study co-mapped and identified distinguishing molecular features of squamous cell carcinomas (SCCs) from five sites associated with smokin

Pan-cancer Alterations of the MYC Oncogene and Its Proximal Network across the Cancer Genome Atlas

Although theMYConcogene has been implicated incancer, a systematic assessment of alterations ofMYC, related transcription factors, and co-regulatoryproteins, forming the proximal MYC network (PMN),across human cancers is lacking. Using computa-tional approaches, we define genomic and proteo-mic features associated with MYC and the PMNacross the 33 cancers of The Cancer Genome Atlas.Pan-cancer, 28% of all samples had at least one ofthe MYC paralogs amplified. In contrast, the MYCantagonists MGA and MNT were the most frequentlymutated or deleted members, proposing a roleas tumor suppressors.MYCalterations were mutu-ally exclusive withPIK3CA,PTEN,APC,orBRAFalterations, suggesting that MYC is a distinct onco-genic driver. Expression analysis revealed MYC-associated pathways in tumor subtypes, such asimmune response and growth factor signaling; chro-matin, translation, and DNA replication/repair wereconserved pan-cancer. This analysis reveals insightsinto MYC biology and is a reference for biomarkersand therapeutics for cancers with alterations ofMYC or the PMN

Spatial Organization and Molecular Correlation of Tumor-Infiltrating Lymphocytes Using Deep Learning on Pathology Images

Beyond sample curation and basic pathologic characterization, the digitized H&E-stained images of TCGA samples remain underutilized. To highlight this resource, we present mappings of tumorinfiltrating lymphocytes (TILs) based on H&E images from 13 TCGA tumor types. These TIL maps are derived through computational staining using a convolutional neural network trained to classify patches of images. Affinity propagation revealed local spatial structure in TIL patterns and correlation with overall survival. TIL map structural patterns were grouped using standard histopathological parameters. These patterns are enriched in particular T cell subpopulations derived from molecular measures. TIL densities and spatial structure were differentially enriched among tumor types, immune subtypes, and tumor molecular subtypes, implying that spatial infiltrate state could reflect particular tumor cell aberration states. Obtaining spatial lymphocytic patterns linked to the rich genomic characterization of TCGA samples demonstrates one use for the TCGA image archives with insights into the tumor-immune microenvironment

Mapping the cellular and molecular heterogeneity of normal and malignant breast tissues and cultured cell lines

Introduction: Normal and neoplastic breast tissues are comprised of heterogeneous populations of epithelial cells exhibiting various degrees of maturation and differentiation. While cultured cell lines have been derived from both normal and malignant tissues, it remains unclear to what extent they retain similar levels of differentiation and heterogeneity as that found within breast tissues. Methods: We used 12 reduction mammoplasty tissues, 15 primary breast cancer tissues, and 20 human breast epithelial cell lines (16 cancer lines, 4 normal lines) to perform flow cytometry for CD44, CD24, epithelial cell adhesion molecule (EpCAM), and CD49f expression, as well as immunohistochemistry, and in vivo tumor xenograft formation studies to extensively analyze the molecular and cellular characteristics of breast epithelial cell lineages. Results: Human breast tissues contain four distinguishable epithelial differentiation states (two luminal phenotypes and two basal phenotypes) that differ on the basis of CD24, EpCAM and CD49f expression. Primary human breast cancer tissues also contain these four cellular states, but in altered proportions compared to normal tissues. In contrast, cultured cancer cell lines are enriched for rare basal and mesenchymal epithelial phenotypes, which are normally present in small numbers within human tissues. Similarly, cultured normal human mammary epithelial cell lines are enriched for rare basal and mesenchymal phenotypes that represent a minor fraction of cells within reduction mammoplasty tissues. Furthermore, although normal human mammary epithelial cell lines exhibit features of bi-potent progenitor cells they are unable to differentiate into mature luminal breast epithelial cells under standard culture conditions. Conclusions: As a group breast cancer cell lines represent the heterogeneity of human breast tumors, but individually they exhibit increased lineage-restricted profiles that fall short of truly representing the intratumoral heterogeneity of individual breast tumors. Additionally, normal human mammary epithelial cell lines fail to retain much of the cellular diversity found in human breast tissues and are enriched for differentiation states that are a minority in breast tissues, although they do exhibit features of bi-potent basal progenitor cells. These findings suggest that collections of cell lines representing multiple cell types can be used to model the cellular heterogeneity of tissues

Dynamic Reprogramming of the Kinome in Response to Targeted MEK Inhibition in Triple-Negative Breast Cancer

Kinase inhibitors have limited success in cancer treatment because tumors circumvent their action. Using a quantitative proteomics approach, we assessed kinome activity in response to MEK inhibition in triple negative breast cancer (TNBC) cells and genetically engineered mice (GEMMs). MEK inhibition caused acute ERK activity loss, resulting in rapid c-Myc degradation that induced expression and activation of several receptor tyrosine kinases (RTKs). RNAi knockdown of ERK or c-Myc mimicked RTK induction by MEK inhibitors, whereas prevention of proteasomal c-Myc degradation blocked kinome reprogramming. MEK inhibitor-induced RTK stimulation overcame MEK2 but not MEK1 inhibition, reactivating ERK and producing drug resistance. The C3Tag GEMM for TNBC similarly induced RTKs in response to MEK inhibition. The inhibitor-induced RTK profile suggested a kinase inhibitor combination therapy that produced GEMM tumor apoptosis and regression where single agents were ineffective. This approach defines mechanisms of drug resistance, allowing rational design of combination therapies for cancer

Integrated Molecular Characterization of Uterine Carcinosarcoma

SummaryWe performed genomic, epigenomic, transcriptomic, and proteomic characterizations of uterine carcinosarcomas (UCSs). Cohort samples had extensive copy-number alterations and highly recurrent somatic mutations. Frequent mutations were found in TP53, PTEN, PIK3CA, PPP2R1A, FBXW7, and KRAS, similar to endometrioid and serous uterine carcinomas. Transcriptome sequencing identified a strong epithelial-to-mesenchymal transition (EMT) gene signature in a subset of cases that was attributable to epigenetic alterations at microRNA promoters. The range of EMT scores in UCS was the largest among all tumor types studied via The Cancer Genome Atlas. UCSs shared proteomic features with gynecologic carcinomas and sarcomas with intermediate EMT features. Multiple somatic mutations and copy-number alterations in genes that are therapeutic targets were identified

Comprehensive Molecular Characterization of Pheochromocytoma and Paraganglioma

SummaryWe report a comprehensive molecular characterization of pheochromocytomas and paragangliomas (PCCs/PGLs), a rare tumor type. Multi-platform integration revealed that PCCs/PGLs are driven by diverse alterations affecting multiple genes and pathways. Pathogenic germline mutations occurred in eight PCC/PGL susceptibility genes. We identified CSDE1 as a somatically mutated driver gene, complementing four known drivers (HRAS, RET, EPAS1, and NF1). We also discovered fusion genes in PCCs/PGLs, involving MAML3, BRAF, NGFR, and NF1. Integrated analysis classified PCCs/PGLs into four molecularly defined groups: a kinase signaling subtype, a pseudohypoxia subtype, a Wnt-altered subtype, driven by MAML3 and CSDE1, and a cortical admixture subtype. Correlates of metastatic PCCs/PGLs included the MAML3 fusion gene. This integrated molecular characterization provides a comprehensive foundation for developing PCC/PGL precision medicine

Integrative Genomic Analysis of Cholangiocarcinoma Identifies Distinct IDH -Mutant Molecular Profiles

Cholangiocarcinoma (CCA) is an aggressive malignancy of the bile ducts, with poor prognosis and limited treatment options. Here, we describe the integrated analysis of somatic mutations, RNA expression, copy number, and DNA methylation by The Cancer Genome Atlas of a set of predominantly intrahepatic CCA cases and propose a molecular classification scheme. We identified an IDH mutant-enriched subtype with distinct molecular features including low expression of chromatin modifiers, elevated expression of mitochondrial genes, and increased mitochondrial DNA copy number. Leveraging the multi-platform data, we observed that ARID1A exhibited DNA hypermethylation and decreased expression in the IDH mutant subtype. More broadly, we found that IDH mutations are associated with an expanded histological spectrum of liver tumors with molecular features that stratify with CCA. Our studies reveal insights into the molecular pathogenesis and heterogeneity of cholangiocarcinoma and provide classification information of potential therapeutic significance